ABSTRACT
We introduce new differentiation matrices based on the pseudospectral collocation method. Monic Chebyshev polynomials (MCPs) were used as trial functions in differentiation matrices (D-matrices). Those matrices have been used to approximate the solutions of higher-order ordinary differential equations (H-ODEs). Two techniques will be used in this work. The first technique is a direct approximation of the H-ODE. While the second technique depends on transforming the H-ODE into a system of lower order ODEs. We discuss the error analysis of these D-matrices in-depth. Also, the approximation and truncation error convergence have been presented to improve the error analysis. Some numerical test functions and examples are illustrated to show the constructed D-matrices' efficiency and accuracy.